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In this paper we apply the formalism of Accelerated Quantum Dynamics (AQD) to the radiative 
stopping of highly relativistic electrons in ice. We compute the acceleration profile of the electron 
along with its lifetime to decay into a muon. The Planckian spectrum of the emitted muon along 
with the its generalized displacement law are presented and used to quantify the muons properties. 

The results predict the acceleration-induced decay of electrons at IceCube energies. The signal 
of electron decay at IceCube manifests itself as an excess of track topologies in an energy window 
accessible experimentally. This setting has the potential to probe the Unruh effect as well investigate 
the flavor content of cosmic ray neutrinos. 
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I. INTRODUCTION 

The dynamics of quantum fields propagating and interacting in classical general relativistic backgrounds is described 
by quantum field theory in curved spacetime. Within its formalism is the natural extension of Minkowskian quantum 
field theory to curved spaces as well as the three canonical particle production mechanisms of the Parker [1], Hawking 
[2], and Unruh [3] effects. With the Unruh effect, there is a production of thermalized particles from an apparent 
horizon created by moving from an inertial into an accelerated reference frame. One can then use this thermal bath 
produced by the acceleration to induce particle transitions [4-9]. This transition rate of particles under acceleration 
is equivalent to the interaction rate of particles in a thermal bath at the accelerated temperature t a = AQD 
[9] is a formalism capable of analyzing a wide variety of these acceleration-induced processes and the observables 
which characterize them. Here we apply the formalism to the enormous deceleration felt by a charged particle, in the 
radiative regime, stopping in ice. Our analysis will focus specifically on the accelerated transition of electrons back 
into muons. We will show there exists an energy scale above which the electrons lifetime is smaller than the time it 
takes to exit the radiative regime and thus has a significant probability for the decay to occur. In the event of an 
electron decay, the energy of the emitted muon, as measured in the proper frame of the electron, has a generalized 
Planck distribution and displacement law that predicts the peak energy to be directly proportional to the accelerated 
temperature. When the muon is boosted back to the lab frame, there exists a large parameter space of energy 
which can be used to investigate the role of acceleration in the muon emission. In particular, there are kinematically 
forbidden regions where the energy of acceleration E a = can be investigated. Moreover, there exists the possibility 
of probing the Unruh effect [3] along with its similarities and differences with the radiative processes of the standard 
model [10]. The analysis is applied specifically to the energies probed by the IceCube detector. The signals of electron 
decay are translated into the signal topologies which characterize the data at IceCube. Evidence for electron decay 
will manifest itself via an increase in the number of track topologies in the energy window where the decay is expected 
to occur. Measurement of the shower to track ratio both in and out of the electron decay window also provides an 
alternative method to directly measure the incoming neutrino flavor content. 

In this paper, Sec. II outlines the use of radiative energy loss in determining the deceleration profile of the particle. 
The relevant kinematic quantities which characterize the motion are also presented. Section III examines electrons and 
muons radiatively stopping in various materials and the energy scales which characterize their motion, i.e. radiation 
length, critical energy, accelerated temperature, etc. A particular emphasis is placed on an ice medium. In Sec. IV we 
compute the electron lifetime and compare it to the other time scales which characterize the motion. We show there 
exists an energy regime where the acceleration can be treated as time-independent and the electron has a statistically 
likely probability of decaying. Section V computes the muon spectrum and determines the most probable energy via 
application of the displacement law. We then boost the muon back into the lab frame and present the parameter 
space of the correlated electron and muon energies. Discussions on the Unruh effect and conservation of energy are 
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included. Section VI outlines the signal topologies at IceCube and how they can be used to measure the presence of 
electron decay in the form of an excess of tracks in the detector. In Sec. VII we summarize the main conclusions of 
the paper. The natural units of H = c = ks = 1 are implemented throughout. 


II. RADIATIVE DECELERATION PROFILE 

In this section we analyze the acceleration profile of a charged particle stopping in a dielectric medium. For ultra 
relativistic electrons, in the radiative regime, the energy loss is dominated by the emission of bremsstrahlung [11]. The 
ratio of radiative energy loss to the energy loss of collisions grows linearly with the electrons energy E. For energies 
above the critical energy, where the two energy loss mechanisms are equal, the radiative energy loss will dominate. 
For electrons in ice, the critical energy is E c ~ 79 MeV. In this manuscript we are considering electron energies on 
the order of E ~ 100 TeV. Thus for our present concern the bremsstrahlung emission dominates over collision events 
by 6 orders of magnitude. This ensures a smooth energy loss and subsequent deceleration. The background medium 
may also induce a transition in an accelerated system, in this case the electron, by coupling to the lattice vibrations 
of the bulk crystalline matter. However, it has been shown [12] that an accelerated detector moving relativistically in 
a background thermal bath will undergo transitions dominated by the accelerated temperature. This implies that any 
coupling between the accelerated electron and vibrations within the ice, i.e. phonons, will be suppressed, and thus 
negligible, in the relativistic regime. The simultaneous emission of Larmor radiation, in addition to Unruh radiation, 
has also been investigated in the cases of electrons in vacuum as well as a plasma [13,14]. The angular distributions 
of the emitted radiation have certain regimes, namely forward and backward proper frame emission directions, where 
they do not overlap and interfere. This ensures the proposed process of electron decay will not be suppressed by 
interference with the emitted bremsstrahlung. Finally, the use of the interactions with bulk matter as a mechanism 
of deceleration to induce particle emission has been studied using laser pulses fired into a dielectric medium [15]. 
In this system the index of refraction, which is characterized by the laser frequency, material properties, and their 
interaction, is responsible for the deceleration. What is of particular importance is that the emitted photons were 
measured experimentally and thus provided confirmation of the theory. In regards to the present manuscript, this 
implies that the use of the interaction of the electron with the background ice, as a method of deceleration, will not 
spoil the acceleration-induced emission process. The bremsstrahlung itself already has a description from the point 
of view of the Unruh effect [16] and its presence alone may indeed a priori imply the validity the analysis. Thus, to 
develop the acceleration profile, we recall the standard energy loss [ 10 ] of charged particles in the radiative regime is 
given by 


dE E 

dx xq 

=>E{x) = E 0 e~ x/x °. (1) 

The radiation length xo determines the distance traveled until the particle loses e -1 of its energy. Here we will only 
consider ultra-relativistic particles with total energy is given by E = 7717 . In this regime, the proper velocity u can be 
approximated via u ~ 7 . The above expression then gives the proper velocity as a function lab frame distance. Thus, 


u(x) = u 0 e~ x/x °. 


( 2 ) 


From this we can separate variables and determine 
Hence 


the lab frame distance traversed as a function of proper time. 


dx 

dr 



=>• x(r) 


■u 0 e 


— x/xo 


J uodr 

x 0 In (t/t 0 + 1 ). 


( 3 ) 


We have defined to = xo/uq and fixed the penetration into the medium at x(t = 0) = 0. The proper velocity and 
acceleration as a function of proper time can then be easily obtained by differentiation. Thus 
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it(r) 

a(r) 


u 0 


(^ + !) 

ao 


(i + !) 2 ' 


(4) 


We have dropped the sign associated with the deceleration since we are only concerned with the magnitude. Also, 
we assume the electron propagates semi-classically along the above derived trajectories, thus giving a well defined 
world line. The acceleration scale is set by the quantity ao = Uq/xq which is both a property of the material as well as 
the initial energy. Recalling for ultra-relativistic velocities uo ~ 7 = Eo/m, we list the relevant kinematic quantities 
which set the scale of our system as follows: 


a 0 


u 0 


to 


xq m 2 

Eo 

m 

Xo m 


(5) 


The critical energy E c of a material is the scale below which the radiative processes no longer dominate. As such, 
we can also define the low acceleration cutoff a c = Xa ^ n i below which we will no longer be able to use the radiative 
stopping power as given in Eqn. (1). We can also define a critical time r c which characterizes the amount of proper 
time it takes for a particles energy to reduce to the critical energy and exit the radiative regime. Inverting the proper 
acceleration in Eqn. (4) and evaluating the proper time at the critical acceleration we have 


t c 



mx 0 


Eo 
mx 0 


Eq 

E c 


1 


( 6 ) 


In the last line we assumed Eq E c . Finally, the radiative energy loss of electrons will dominate up until the energy 
scale where the LPM effect begins to spoil the deceleration [17]. This suppression of photons with wavelengths of the 
order of the distance between scattering sites begins to diminish the energy loss near the LPM energy scale Elpm- 
The above deceleration profile will therefore only be valid for incident electrons with energy Eq < Eppm- Under the 
LPM threshold, the relevant parameters that characterize the system can be written in terms of the incident particle 
energy, incident particle mass, the radiation length, and critical energy of the decelerating material. We will now 
analyze the relevant kinematic quantities specifically for electrons and muons in various materials with a particular 
emphasis on ice. 


III. ENERGY SCALES OF ELECTRONS IN ICE 

In the previous section we outlined the relevant parameters which characterize our decelerated motion. Table 
I contains the relevant kinematic parameters for various materials. Moreover, the focus of this paper will be the 
excitation of electrons back into muons and, as we shall see in the next section, the energy scale of the acceleration 
is best written in terms of the muon mass. Defining the dimensionless acceleration a 0 = ao/m M , and likewise for the 
critical acceleration, we can look at the relevant energy scales for various materials. To compute the accelerations, we 
assume an incident electron energy of Eq = 100 TeV. 
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Material 

xo [cm] 

E c [MeV] 

Elpm [TeV] 

«c [xlO~ 10 ] 

5o 

Uranium 

0.317 

6.65 

2.39 

1 

22700 

Lead 

0.561 

7.43 

4.32 

.708 

12830 

Iron 

1.757 

21.7 

13.5 

1.93 

4096 

Ice 

39.31 

78.6 

303 

1.13 

183.1 



TABLE I: The radiation length, critical energy, LPM energy, critical acceleration, and initial acceleration for 100 TeV electrons 
in various materials. 

Focusing specifically on ice, Table II contains the dimensionless acceleration scale for different incident electron 
energies. For completeness, we include the accelerated temperature t a which also characterizes the overall energy 
scale of the decay process. The dimensionless acceleration enables these acceleration-induced processes to occur and 
regimes where do ^ 1 are where these effects are strongest. Moreover, in the interest of measuring these effects, ice is 
the easiest since it is transparent to light and thus facilitates measuring the energy of the system via bremsstrahlung. 
It is worth noting that for incident electron energies E 0 < 500 TeV the accelerated temperature t a is at, or below, the 
weak scale t a < m w ±. We leave an analysis that includes the weak gauge bosons at higher energy scales for future 
work. 


E 0 [TeV] 

5o 

t a [GeV] 

10 

1.831 

.031 

50 

45.77 

.77 

100 

183.1 

3.08 

300 

1637 

27.5 

500 

4577 

77 



TABLE II: The acceleration scale and temperature for electrons at various energies in ice. 

It is also worth tabulating the same parameters for muons in various materials. Table III looks at the kinematic 
parameters of 100 TeV muons moving through matter in the radiative regime. Here we use the same dimensionless 
acceleration written in terms of the muon mass. It should be noted that muons have a much larger radiation length 
and critical energy than electrons. Moreover the LPM energy scale for muons is a factor ~ 10 10 larger than electrons 
and is therefore negligible [10]. This allows us to analyze certain energy regimes where the electron is in the radiative 
regime and decay is most probable but the final state muon does not decay. This will better enable us to accurately 
analyze these processes. 


Material 

xq [nr] 

E c [GeV] 

d c [x 10 -10 ] 

a 0 [xl0~ 5 ] 

Uranium 

26.69 

128 

1.046 

6.387 

Lead 

48.03 

141 

.706 

3.5497 

Iron 

152.5 

347 

1.346 

1.1178 

Ice 

3120 

1031 

.58074 

.054634 



TABLE III: The radiation length, critical energy, critical acceleration, and initial acceleration for 100 TeV muons in various 
materials. 

Then, for ice, Table IV contains the dimensionless acceleration and accelerated temperature for muons. It should 
be noted that the acceleration scale and temperature are substantially smaller for the muon than for the electron. 
This is the result of having a much larger radiation length and critical energy. Moreover this coincidence helps to 
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ensure that the accelerated decay processes we are looking for are more prominent for electrons than for muons. It 
should also be noted that the accelerated temperature for muons is substantially lower than the weak scale and the 
acceleration scale do 1 . 


E 0 [TeV] 

a 0 [x 10 -7 ] 

t a [KeV] 

100 

5.4634 

.009 

200 

21.8536 

.038 

400 

87.4146 

.147 

800 

349.658 

.588 

1000 

546.341 

.919 



TABLE IV: The acceleration scale and temperature for muons in ice at various energies. 

In this section we tabulated the various parameters that characterize the deceleration of both electrons and muons 
in various materials. A particular emphasis was placed on ice so we can search for these processes at the energy scales 
of IceCube [18]. The transparency of ice to light, as well as the distinct topologies of electrons and muons will help 
characterize the energies as well as the type of event that occurs [19]. We saw that there isn’t sufficient stopping 
power in the ice to bring the muons to the necessary acceleration scale. The electrons however do have sufficient 
deceleration. It is this asymmetry that will help define a possible event. 


IV. THE ELECTRON LIFETIME 

The formalism of AQD [9] is capable of computing a wide variety of observables of acceleration-induced interactions. 
In this section we compute the lifetime of an accelerated electron. This process is essentially the inverse beta decay 
of muons rather than protons. The relevant transition is given by 

e ± T„/l ± + f (1 + l'e. (7) 

Here, we will model all final state particles as being massless since the accelerated temperature of the system is 
much larger than the muon and neutrino rest masses. The neutrinos can be easily approximated as massless particles. 
The acceleration-induced transition rate for arbitrary n-particle final state multiplicity [ 8 , 9] is given by 


T(A E,a,n)=G 2 n 



1 


(An — 2)!! 


n—1 

n 

k—0 


1 + k 2 



1 

e 2-KAE/\a\ _ 1 ’ 


( 8 ) 


We note this formalism does include time dependent accelerations. Now, to model the electron decay we note 
there will be three final state particles, i.e. n — 3 for the muon and two neutrinos. Moreover, since the electron 
and neutrino masses are much smaller than the muon, the transition energy will be the muon mass, i.e. A E = 
+ m v + m Ve — m e ~ m /i . The electron decay rate is then given by 


r (n) = r 2 ™l 1 + 5d 2 + 4d 4 
38407T 5 e 27 I 7l“l — 1 


(9) 


Note, we have written the decay rate in terms of the dimensionless acceleration a = a/m For transitions between 
the electron and muon the energy scale is set by the muon mass. We can fix the coupling G 2 via taking the inertial limit 
of the acceleration-induced muon decay and use the detailed balance to import the coupling for the electron excitation 
[ 8 ]. The acceleration-dependent muon decay rate can be obtained from Eqn. (9) simply by letting m M —> —to^, i.e. 
we transition down in energy. The daughter electron and neutrinos can be considered massless which again yields 
n = 3. This technique is equivalent to using the detailed balance of the two process at thermal equilibrium. Also note 
that this requires |a| —>• — |a|. The acceleration-dependent muon decay rate [4,8] is then given by 
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r^( 5 ) 


„2 1 + 5 a 2 + 4a 4 

38407T 5 1 — e _27r /l a l 


The inertial muon decay rate [10] is well known to be 


( 10 ) 


\ _ G >1 
^ ~ 1927T 3 • 


( 11 ) 


We can then fix the coupling by taking the inertial limit of the acceleration-dependent muon decay rate and matching 
it to the inertial rate. Thus 


lim T Al (d) = A p 


a—± 0 


35 


=> G = —7r J Gi. 


( 12 ) 


Now that we have fixed the coupling, we have the acceleration-induced electron decay rate is given by 


r e (a) = 


G 1 + 5 a 2 + 4 a 4 

1927T 3 e 27r/|a.| _ l 


(13) 


If we recall the canonical muon lifetime is given by = 2.184 /ns, we can determine the lifetime by merely 

reciprocating the decay rate. Hence 


T e (a) 


e 2 ^/a _ x 
Tfl 1 + 5a 2 + 4a 4 ' 


(14) 


To better understand the decay process we must examine the various time scales of the system. We know the 
critical time r c from Eqn. ( 6 ) determines how long the electron would take to exit the radiative regime under the 
assumption of an initial energy much larger than the critical energy. For electrons in ice the critical time is given by 
r c = m e?o — g 52 x 10 ~ 12 s. Furthermore for lifetimes smaller than the time scale tq from Eqn. (5), we can assume 

the acceleration to be constant. For our system the time scale is given by tq = m ^ X( > = 6 ' 7x J ,° — s. Here we are 
measuring Eq in units of TeV. Finally, we must examine the lifetime r e in the limit oflrigh acceleration, i.e. a ?n M . 
This will enable us to determine if the electron will decay prior to exiting the radiative regime, i.e. r e < r c . Recalling 
the form of the time-dependent acceleration from Eqn. (4) we can find the deceleration time 74 , i.e. the maximum 
proper time that can elapse while still allowing the electron lifetime to be smaller than the critical time. Hence 


t c > 

> 

> 

=> T d < 


e 2^/a _ X 

Tm 1 + 5d 2 + 4a 4 
7 r 1 
T/i 2 ^ 

7 Tm^ijdJro + l ) 10 


ro 


t c 2 Eq° 1 
r M 7 r m 40 m^x g 


1/10 


- t 0 . 


(15) 


The above deceleration time evaluates to 74 < 2.5 x 10 -1 ' — 6,7x i ) 4 — s. For lifetimes less than the deceleration 
time, the electron will have the maximum probability to decay before it exits the radiative regime. Now, to analyze 
the lifetime of the electron we will look at the case of constant acceleration. As such we focus on the acceleration 
scale ao■ This will provide a first estimate for the electron lifetime and enable us to analyze the system using the 
initial electron energy E 0 . We note m M ?n 2 Xo = 55 TeV . Then, we find a 0 = ram ° m = -gf- Note we are writing the 
initial energy in units of TeV. Our electron lifetime is then given by 
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Te(Eo) 


e 2n55/E% _ 1 


1 + 5 




(16) 


In order to estimate if the electron decay has enough time to take place, we require the electron lifetime to be 
shorter than the deceleration time, i.e. the electron is statistically likely to decay before it leaves the radiative regime. 
Moreover, if the lifetime of the electron is smaller than the time scale tq than we can disregard the time-dependence 
of the acceleration in the analysis and use the acceleration scale clq to characterize the system. In Fig. 1 we plot 
all relevant time parameters which characterize the decelerating electron. It should be recalled that the lifetime of 
a decaying particle defines the characteristic time for an ensemble to lose - of the initial population. The statistical 
nature of the decay should be kept in mind as we estimate the probability of decay based on the relevant parameters 
of the system. 



6[TeV] 


FIG. 1: Comparison of electron lifetime, critical time, deceleration time, and time scale as a function of the initial electron 
energy. 


The shaded region of Fig. 1 indicates the regions where the lifetime is smaller than the deceleration time and we 
expect an electron decay to occur. Moreover, by inspection we see the energy scale at which the lifetime becomes 
smaller than the deceleration time is also very close to the energy scale time below which we may accurately approxi¬ 
mate the acceleration as constant. We can determine this threshold energy E t for the lifetime to be shorter than the 
scale time tq . Using a time-independent acceleration we find the threshold energy in the same way we computed the 
deceleration time in Eqn. (15). Thus, in the limit of high acceleration we have 




e 27r / a - 1 

TO 

> 

7)1 1 + 5a 2 + 4a 4 



7r 1 


> 

T/i 2 d? 

x 0 m e 

> 

_ 

E 0 

T/i 2 Eq° 

=+ Eo 

> 

7T e: n 4 1 

7 >2 r7 V neX ° 


(17) 


Then, for our electron in ice system, we find that for Eo > 111.1 TeV the acceleration is effectively constant while 
simultaneously satisfying the condition that the electron lifetime is less than the deceleration time. It is above this 
energy scale we expect to see an electron decay event. We can also calculate the final energy of the electron at the 
time of the decay using the proper velocity in Eqn. (4). Figure 2 contains the electron incident energies Eo , the 
energy at the time of decay Ef , and the change in energy Ed = Eq — Ef. At IceCube, the change in energy is, up to 
detector efficiency and effective volume, the measured energy that is deposited into the detector. 


















FIG. 2: Comparisons of the initial energy, final energy, and deposited energy in ice. 


It is interesting to note that above the threshold energy E t ~ 111.1 TeV the electrons do not lose any appreciable 
fraction of energy prior to converting to muons. However, at the ~ 100 TeV scale there will still be a measurable, 
and quite large, amount of energy deposited into the detector. We also note the condition that the time-dependent 
formalism used in the analysis remains valid is that j/a 2 -C 1, where j is the jerk [9]. We know the time-dependence 
of the acceleration becomes irrelevant for energies greater than ~ 111.1 TeV. However we must still verify the time- 
dependent formalism is valid to begin with. By taking the derivative of the proper acceleration we have j = , 
with jo = oo/tq. Then the constraint implies 


J 


< 1 


E 0 > 


2 m e 


[1 - -Z-] 

L TnUn J 


(18) 


Pertaining to this analysis we recall that our lifetimes are much smaller than the scale time To- As such, we require 
Eq 1.035 MeV, which is most assuredly satisfied for all energies under consideration. Thus the condition necessary 
to ensure a valid time-dependent formalism has been verified. We close this section by commenting on any effect 
due to the sharp turn-on of the deceleration when the electron is produced [20]. For an acceleration turn-on time 5, 
energy gap A E, and lifetime r the conditions for a thermal spectrum are r S and r 1/A E. The first condition 
ensures that any transients due to the turn-on will be damped away and the second condition ensures that acceleration 
occurs long enough for the system to thermalize. In our analysis the turn-on time is given by the time it takes for 
the electron to traverse two scattering sites of ice. For ice, the lattice constant is l s ~ 2 x 10“ 10 m and the turn-on 
time is given by <5 = £ s /it 0 = f s T 0 . Note we have defined the fraction f s = £ s /xo ~ 5 x 10 -10 . As a function of the 
initial energy we have / s tq ps 3x ^ — s. Within the energy window analyzed here, Eq ~ 100 TeV, the turn-on time is 
given by / s 7o w 3 x 10~ 23 . From Fig. 1 we clearly see that 10~ 18 s > r e > 10“ 22 s and therefore the condition that 
the lifetime is greater than the transient period, r e / s ro, is satisfied within the energy window where the decay is 
expected to occur. To ensure tliermalization, we note the energy gap for the transition is given by the muon mass, 
A E = m^. The condition, r e l/m M , can easily be verified by noting 1 /m M ~6x 10 -24 s. Therefore, based on the 
electron lifetime constraint from Fig. 1, the system has sufficient time to thermalize, i.e. r e l/m M . Thus, for the 
acceleration-induced electron decay at IceCube, the conditions which are necessary to produce a thermal spectrum in 
the presence of a sharp turn-on are satisfied. 


V. THE MUON SPECTRUM 

Having determined that there is indeed an energy regime in which we expect an electron decay to occur, we now 
endeavor to compute the energy of the emitted muon. Although there are also neutrinos emitted in the decay process, 
we exclude them from the analysis on account of their leaving the detector with virtually zero interaction with the ice. 
The muons, however, are charged and will deposit energy, via bremsstrahlung, as they propagate away. Moreover, as 
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we found in the first section, muons at IceCube energies will not have any appreciable deceleration in comparison to 
the electrons. The generalized spectra Af from [9] is given by 


AT{AE, a, n) 


1 G 2 w / A E + Co 

T (27r) 2 \ 7T 


2n—3 ^ n—2 


(4n — 6)!! 


n 

k—0 


1 + fc 2 


A E + Cj 


1 


g27r(AS+cIj)/|a| _ ^ ' 


(19) 


The parameter Co is the energy of the emitted particle, in this case the muon, as measured in a comoving frame 
instantaneously at rest with the accelerated particle, in this case the electron. Enforcing n = 3 for the electron 
excitation, and imposing the coupling G 2 = ^7r 2 G 2 we have 


Afn(a,Co) = T e (a) 


ttG^Coirn^ + w ) 3 


1 


( —-—V 

\ m^+uj J 


15367T 3 g27r(m (1 +il>)/|a| _ ^ 

Then writing everything in terms of the initial energy Eq we find 


( 20 ) 


AT/j ( E 0 , Co) 


■XEo)- 


35 G 2 




-N 3 i + ( M V 

w ) 1 + \m'ixoim^+ij) ) 


15367T 3 


e 2'irmlx 0 (m IJ .+u)/E^ _ ^ ' 


( 21 ) 


Figure 3 contains spectra for various initial electron energies normalized to unity over Co via N = j ■ Recalling 
the threshold energy is E t = 111.1 TeV, we will restrict our analysis to energies above the threshold so as to ensure the 
appropriate time-independent analysis. The Planckian spectra reflect the thermal nature of the accelerated reference 
frame. 



FIG. 3: The spectrum of the emitted muon for various initial electron energies. 

The notoriously low statistics associated with neutrino interactions necessitates we examine the most probable 
energy emitted rather than the distribution. This is accomplished via determination of the peak energy of the emitted 
muon spectra using the generalized displacement law [9] 


xe 


e x - 1 


1 - 


2ttAE 

\a\x 


+ (2 n -3)-2(^) Y, 


2 n —2 




= 0 . 


( 22 ) 


The displacement coefficient x is computed numerically and determines the peak muon energy via the expression 


Co = x- -A E. 

2 tt 


( 23 ) 
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Here we find the inherently quantum mechanical energy of acceleration E a = The forthcoming analysis shows 
there exists the potential to experimentally probe the nature of this energy. Now, application of the displacement law 

to our system necessitates enforcing the condition that n = 3, A E = m^, and a = ao = m a Xg ■ As such, we obtain 


xe x 
e x — 1 


+ 3~ 



0 . 


(24) 


Numerically solving for x , for any incident electron energy, enables us to then determine the most probable muon 
energy as measured in the proper frame instantaneously at rest with the accelerated electron. It is worth mentioning 
that the muons energy comes from the accelerated temperature. The peak muon energy wo is given by 


w 0 




, El 

2-7t55 



(25) 


In Table. VI we have tabulated the displacement coefficient for various initial energies. We see that to within a 
percent we have x ~ 1.8. It is also interesting to note that the proper energy of the emitted muon grows as the square 
of the initial electron energy albeit weighted by a small prefactor. 


E 0 [TeV] 

X 

wo [GeV] 

111.1 

1.818 

6.758 

200 

1.800 

21.92 

303 

1.796 

50.33 



TABLE V: The displacement parameter and associated peak muon energy for each incident electron energy. 

We can finalize this analysis by noting that emitted muon will also be boosted due to any residual velocity from the 
incident electron. We can now use the electron lifetime to compute the electrons proper velocity and thus the Lorentz 
gamma. We also compute the Lorentz gamma of the muon as measured in the electrons frame. This will enable us to 
determine the energy of the muon as measured in the lab frame. First, we know the muons energy, in the electrons 
frame, is ujq = From the muons Lorentz gamma we can also invert it to get the corresponding proper velocity. 

Hence 


7m = 




h 


Wo 

% 



(26) 


For the electron we know at the time of emission q e = and we can also solve for the electrons velocity f3 e = 
y/l — l/ 7 g. Then using the composition of Lorentz gammas via the relativistic velocity addition formula [11] 7 ^ = 
[1 + /3 e /3 p cos (0 M )]TeT^ we have the muon energy, as measured in the lab frame, is given by w = The angle 9 M 

is measured in the proper frame of the electron and determines the angle between the muons velocity relative to the 
electrons velocity. For both forward and backwards emitted muons we have 


w± 


rrif, 7 ± 

± PeP/ihelv 


Wo 


El 

m e 


1 ± 


1 - 1 ^ 


Ef 


1/2 



(27) 


Note, the plus or minus corresponds to muon velocity being parallel or anti-parallel to the electrons velocity 
respectively. These limits bound the muons energy when measured in the lab frame. Figure 5 contains a plot which 



11 


outlines all emitted muon energies as a function of the electrons incident energy using the approximate displacement 
constant x = 1.8 for all energies. The shaded region to the right of the threshold energy E t partition comprises 
the domain of muon energies predicted by AQD. The forward and back scattered limits are denoted by w+ and ir¬ 
respectively. Regions above the initial electron energy are kinematically forbidden in Minkowski space and could 
provide the opportunity to investigate the nature of the energy of acceleration E a . Below the initial electron energy 
both the standard model and AQD can describe the signal. Here an analysis comparing the energies of muons emitted 
by other processes would be needed to differentiate the signal. This region may also be used to compare the emission 
rates, and energies, of muons from the Unruh effect and from standard model radiative processes. 



FIG. 4: Comparison of the initial electron energy and final muon energy as measured in the lab frame. 


The parameter space of muon energies which kinematically conserve energy can be used to constrain the muon 
emission angle. For ultra relativistic electrons and muons /3 e /3„ ~ 1, and the final lab frame muon energy w M = Ef 
we have 1 = ^-[1 + cos Inverting the cosine then yields the muon emission angle. Hence 


2m e 


w 0 


1/2 


2 m e 


m. 


( E o 
2-7r55 


— 1 


1/2 


(28) 


Thus for conservation of energy to be imposed, we require the emitted muon to be back scattered in the reference 
frame of the electron. This implies the muon is emitted away from the apparent horizon of the accelerated reference 
frame. After the muon has been emitted it will lose energy via a combination of ionizing and radiative processes [10]. 
This energy loss can be parametrized as 


doj 
dx 
=t- w(x) 


—a — /3w 

uj oe ~ x/xo + ^(e~ x/xo - 1). 


(29) 


For muons in ice we have xq = 1//3 = 3120 m and a/(3 = .4762 TeV. With these parameters, along with the 
monotonically increasing energy deposited by the muon w^ = Wo — w(x), we can determine the initial energy of the 
muon to be 


w 0 


Ud _ a 

l_ e -x/x 0 p- 


(30) 
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The understanding of the energetics of both the electron and muon, in particular their deposited energies, enables 
one to confirm this acceleration-induced process. What is interesting to note is that there exists a rather appreciable 
parameter space where the muon has more energy than the initial electron. To conserve energy, a constraint is made 
on the muon emission angle which implies the muon is emitted away from the horizon. Moreover, any apparent 
violation of conservation of energy could be used to probe the energy of acceleration E a = The correlations 

between the emitted muon energy and the initial electron energy, along with a comparison with background muon 
and scattering emission rates, has the potential to distinguish the nature of the muon emission and thereby probe 
the existence of the Unruh effect. We must also mention the energy loss via bremsstrahlung will be suppressed at 
sufficiently high energies Elpm > 303 TeV due to the LPM effect [17]. Above this energy threshold, the electrons 
energy loss is no longer sufficient to produce the required acceleration and the electron decay probability goes to zero. 


VI. ELECTRON DECAY SIGNAL AT ICECUBE 


When a neutrino enters the IceCube detector it may interact with the ice via charged current or neutral current 
interactions. In neutral current interactions, a Z boson is exchanged between the incident neutrino and a nucleon 
in the ice. The incident neutrino scatters off the nucleon and and deposits energy in the form of a hadronic shower 
produced by the recoil of the nucleon. These signals have spherical topologies in the IceCube detector. In charged 
current interactions, a boson is exchanged between the incident neutrino and a nucleon in the ice. The neutrino 
produces the associated lepton in the final state along with a hadronic shower due to the nucleon recoil. The lepton 
produced then deposits its energy via a combination of radiative and ionization processes depending on the flavor. 
For electrons, the radiative energy loss is of such strength that the electron stops in the detector depositing an 
approximately spherical distribution of energy. For taus, the lifetime is so short that the decay and subsequent 
daughter products are all confined to the detector and deposits a spherical distribution of energy as well. The analysis 
here is done below the energy scale necessary to resolve the so called double bang signals of the creation and subsequent 
decay of the tau. These spherical signals are labeled as showers. For muons, the lifetime is of sufficient length that 
the muon leaves the detector depositing a track of energy deposited along the way. 

The signals are classified by their event topologies; either showers or tracks. For any incoming neutrino, one third 
of the event signals will be due to neutral current interactions and the remaining two thirds of the event signals will 
be due to charged current interactions. For tau neutrinos, this implies all signals will have shower topologies. For 
muon neutrinos, one third of the events will have shower topologies while the remaining two thirds will have track 
topologies. For electron neutrinos one third of the events will have shower topologies. The remaining two thirds of 
the interactions which produce the associated electron will have shower topologies if the electron does not decay or 
track topologies if the electron decays into a muon. If we consider N total incoming neutrinos species with the relative 
fractions f e , f tl , and f T such that f e + / M + f T = 1, we can classify the signal via the shower to track ratio |. In the 
absence of electron decay, the total showers s, tracks t, and the shower to track ratio cr = | is given by 


s 

t 

a 


f e N + -f^N + f T N 


\f, N 

, 3 (fe + /r) 

2 L u 


(31) 


To classify the event topologies in the presence of electron decay we denote each variable with a tilde, i.e. s, t, and 
<7 = |. Then for electron decay we have 


s 

t 

a 


lfeN+^UN + f T N 
ifeN+^UN 


(32) 


Finally, to probe the signal of electron decay, all one needs to analyze is the shower to track ratio at energies inside 
and outside of the electron decay window E t < E < Eppm- A difference in the measured ratios would warrant a 
more detailed search and, provided the flavour content is not energy-dependent, could be explained by electron decay. 
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One could also use the measured ratios to pin down the incoming flavor content. Note Eqn.(31), Eqn. (32), and the 
unitary condition enable us to solve for the incoming flavors. Hence 


fe 

fn 

fr 


3 (<r — a) 

2 (1 + < t ) (1 + ct ) 
3 

2 (1 + ct ) 

2a - 1 
2 (1 + <t) 


(33) 


As an example we note that for incoming neutrino flavors of equal probability, i.e. / e = = / T = |, one would 

measure a = ^ and a = Another way to see this is simply as an excess of tracks with a sharp turn on at 
the threshold energy and a bit of a slower cutoff near the LPM energy. This analysis assumes there is no energy 
dependence in the relative fractions of the incoming neutrinos. 

We close this section by commenting on the inclusion of fermions in the analysis rather than using scalar fields. 
By comparing the polynomials of multiplicity for an n = 2 fermionic description from [2] and an n = 3 description 
used here, we note that the inclusion of fermions yield a higher order polynomial of multiplicity. The effect of this 
in the analysis would be to push the threshold energy Et to a lower energy scale. This would increase the electron 
decay window to encompass a larger number of neutrino signals at IceCube. We also must mention that the energy 
scales are consistent with a valid Fermi theory analysis. The fermionic computation developed in [5-7], adapted to 
electrons and muons, and without the use of an Unruh-DeWitt detector would provide a significantly more accurate 
prediction that satisfies all electroweak t ao <C m w ±, massless do time-independent t e Ct 0 , and perturbative 

approximations G 1 2 3 4 5 6 7 < 1. 


VII. CONCLUSIONS 

In this paper we utilized the radiative deceleration of electrons in ice to induce their decay back into muons. 
The acceleration profile, in the radiative regime, was developed and the relevant energy scales were tabulated for 
electrons and muons in different media including ice. We computed the electron lifetime and showed that a time- 
independent acceleration was sufficient for analysis beyond an accessible threshold energy. We also showed that above 
the energy threshold the electron lifetime was smaller than the critical time to exit the radiative regime and thus had 
an appreciable probability of decaying. The muon spectrum was computed in the electrons proper frame and also 
boosted to the lab frame. Using the radiative energy loss of both the electron and muon, we also computed the energy 
deposited in the ice by both particles. The analysis was carried out at IceCube energies. A description of the signal 
topologies in the presence of electron decay was also included. The resulting analysis shows that IceCube has the 
potential to verify the formalism of AQD, provide experimental evidence for the Unruh effect, investigate the nature 
of the quantum mechanical energy of acceleration E a , and pin down the incoming neutrino flavor content. 
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